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We re-examined the angular dependence of the radiation from the intrinsic Josephson junctions 
in rectangular mesas of Bi2Sr2CaCu2 08+(5, in order to determine if the cavity mode part of the 
radiation arises from waves across the width w or along the length i of the mesas, associated with 
"hot spots" [Wang et a/., Phys. Rev. Lett. 105, 057002 (2010)]. We derived analytic forms for 
the angular dependence expected in both cases for a general cavity mode in which the width of the 
mesa corresponds to an integer multiple of one-half the wavelength of the radiation. Assuming the 
coherent radiation from the ac Josephson current source and the cavity magnetic surface current 
density source combine incoherently, fits to the data of Kadowaki et al. [J. Phys. Soc. Jpn. 79, 
023703 (2010)] on a mesa with mean l/w = 5.17 for both wave directions using two models for the 
incoherent combination were made, which correspond to standing and traveling waves, respectively. 
The results suggest that the combined output from the uniform ac Josephson current source plus a 
cavity wave forming along the rectangle length is equally probable as that of the combined output 
from the uniform ac Josephson current plus a cavity wave across the width. However, for mesas 
in which n^/2w is integral, where n is the index of the rectangular TM^ q mode, it is shown that 
standing cavity mode waves along the length of the mesa do not radiate in the xz plane perpendicular 
to the length of the mesa, suggesting experiments on such mesas could help to resolve the question. 

PACS numbers: 07.57.Hm, 74.50.+r, 85.25.Cp 



I. INTRODUCTION 

Since there is a gap in the range 0.1-10 terahertz (THz) 
in the electromagnetic spectrum for practical coherent 
sources, development of a material as a useful source 
could lead to many practical applications in the detec- 
tion of complex molecules, in national security detection 
issues, space communications, environmental engineer- 
ing, biochemistry analysis, cancer detection, and possi- 
bly genetic engineering [jy] . Recently, the discovery of in- 
tense, continuous, coherent sub-terahertz radiation upon 
the application of a static voltage V across the N intrin- 
sic Josephson junctions in mesas of Bi2Sr2CaCu208+5 
(Bi-2212) has generated a lot of excitement, in the hope 
that such a material may lead to a useful source filling 
the "terahertz gap" [1, [H. As in those papers, it was 
generally found that the fundamental angular frequency 
ujj of the ac Josephson current generated by the applied 
dc voltage V satisfies the ac Josephson relation 
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where e is the magnitude of the electronic charge, % is 
Planck's constant divided by 27r, and the fundamental 
Josephson frequency fj = ujj/{27t). An additional con- 
straint on the radiation occurs when V is adjusted in 
the mesa's current- voltage {I /V) characteristic such that 
the frequency matches that of an electromagnetic cavity 
mode of the rectangular mesa. We assume the rectangu- 



lar mesa has width length i > and height h <^ w^i. 
Writing k and x for the wave vectors and position vectors 
outside the mesa, and k' and x' for the corresponding 
vectors inside the mesa, the angular frequencies ujmp of 
the very thin rectangular cavity TM^ modes are 
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where c is the speed of light in vacuum and rir = ^/e ^ 4.2 
is the index of refraction of Bi-2212 in the frequency range 
of study|2|, |3|. Hence, this additional constraint may be 
written as 
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where mo and po are particular values of m and p. From 
fits to many rectangular mesas [J], it is evident that in 
all studies to date, the frequency /i^o of the TMfg mode 
with mo = l,Po = matches the fundamental frequency 
fj of the radiation. Therefore, most workers assumed 
that the fundamental mode was described by a stand- 
ing wave of wavelength X = 2w across the width of the 
mesa |2|-[l3]- Many theoretical studies were made based 
upon the coupling of the ac Josephson current to this 
particular half-wavelength electromagnetic cavity mode 
[5|-[3- However, with one exception [l6j, these calcu- 
lations were unable to predict any angular dependence 
of the output radiation from the cavity. In the sole ex- 
ception [16], the mesa was modeled as a two-dimensional 
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sheet in either the xy plane or the xz plane, and the out- 
put power was calculated numerically. In that numer- 
ical calculation, it was therefore difficult to obtain the 
full three-dimensional angular dependence of the output 
power for comparison with experiment. Thus, several 
workers employed the magnetic Love equivalence princi- 
ple [19], a variation on Huygens' principle, to replace the 
electric field in the cavity by a magnetic surface current 
density, and then to use that effective magnetic surface 
current density to obtain the three-dimensional far-field 



radiation distribution pjj, [171, |1 



However, in fits to experiment on rectangular mesas, it 
was found that in addition to the radiation from the ex- 
citation of the electromagnetic cavity mode correspond- 
ing to the half- wavelength mode by the non-uniform part 
of the ac Josephson current, the uniform part of the ac 
Josephson current was also found to radiate [4]. This ra- 
diation was treated by the electric Love equivalence prin- 
ciple (l9|, in which the ac Josephson current inside the 
mesa was effecti vely replaced by an electric surface cur- 
rent density [Hf. In fits to the observed output ra- 
diation from a rectangular mesa, it was then found that 
the observed radiation was consistent with a compara- 
ble amount from both the uniform ac Josephson cur- 
rent source and the electromagnetic cavity mode excited 
by the non-uniform part of the ac Josephson current in 
the mesaQ. Direct evidence for the importance of the 
uniform ac Josephson current source was subsequently 
provided by the investigation of cylindrical mesas [2 Qj. 
In that experiment, three disk-shaped cylindrical cavi- 
ties were found to exhibit radiation not only at the fun- 
damental electromagnetic cavity mode for a cylindrical 
cavity, but also at the integral second harmonic of the 
fundamental radiation frequency [20]. Since the modes 
for a cylindrical cavity are not integral multiples of one 
another (in [Til . it was therefore determined that the ob- 
served output radiation was a combination of that arising 
from the excitation of the fundamental cylindrical cavity 
mode excited by the main non-uniform part of the ac 
Josephson current and the radiation from the uniform 
part of the ac Josephson current itself lH, [2o| . The inte- 
grated intensities of the outputs from those two sources 
were found to be comparable to one another [20]. Hence, 
the dual-source mechanism of radiation from both the 
uniform part of the ac Josephson current and from the ex- 
citation of a cavity mode from the appropriate part of the 
non-uniform ac Josephson current was established [20]. 

However, the observation of "hot spots" in rectangular 
Bi-2212 mesas accompanied by the apparent evidence of 
standing waves with \^ 2w along the length of the mesa 
have led to a lot of interest [2l|. Very recently, Wang et 
al. provided evidence that the samples exhibiting such 
hot spots do in fact emit radiation ^]. Although there 
is widespread agreement that the fundamental mode in 
all rectangular mesas studied to date has the frequency 
of the TMf Q mode, it has therefore been an open ques- 



tion as to whether the waves contributing to the radia- 
tion form across the width of the mesa, as expected from 
standard cavity theory 23, [1^, or form with the same 
wavelength \ = 2w along the length of the mesa. In this 
paper, we address that issue in detail, performing de- 
tailed fits to the data of Kadowaki et al. for the angular 
dependence of the radiation in the xz and yz planes nor- 
mal and parallel to the mesa length, respectively ^] , and 
find that both cavity wave directions make equally prob- 
able contributions to the radiation for that mesa with a 
mean l/w = 5.17. However, we show that in the special 
cases in which nl/w = 2g, where g is a natural number 
and n is the index of the TM^ q mode, an excited cavity 
wave formed along the length of the mesa cannot emit 
radiation in the xz plane normal to the mesa length, and 
we thus propose detailed angular dependence studies of 
the radiation from such mesas as a more precise test of 
this question. 

In Sec. II, we present the angular dependence of the 
radiation from a rectangular mesa with the radiation aris- 
ing from two sources: the uniform ac Josephson current 
source, and the radiation from an excited electromagnetic 
cavity mode, in which either standing or traveling waves 
may form either across the width or along the length of 
the mesa. In Sec. HI, we present our fits to the available 
published data on the angular dependence of the radia- 
tion from rectangular cavities. In Sec. IV, we consider 
the special cases of nl/w = 2q. Finally, in Sec. V, we 
summarize and discuss our results. 



II. ANGULAR DEPENDENCE OF THE 
RADIATION 

From previous studies, it was determined that the ra- 
diation from rectangular and cylindrical mesas of Bi-2212 
arose from a dual-source mechanism 0, [20|, in which 
the magnetic field associated with the uniform part of 
the ac Josephson current and the electric field associated 
with the nth ac Josephson current harmonic at fn = nfj 
and the nth cavity mode of the mesa were treated us- 
ing the Love equivalence principles as electric and mag- 
netic surface current sources Jsnix^^t) and Msn{x^ ^t), 
respectively 0, Il3^20, 23, 24]. In rectangular source co- 
ordinates {x'^y'^z') for the directions along the width, 
length, and height of the mesa, respectively, for an elec- 
tric surface current density Jsn{'^') and either a magnetic 
surface current density Msn,w{'x') or Msn,£{x^) for the 
nth mode forming across the width or along the length 
of the mesa, respectively. 
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Msn,e{x') = -^■r]M{z')sm[n{y' - yn)n/w] 

xY,a[y'Ux',y')-x'gAx',y% (6) 

Uix',y') = wSix' + aw/2)e[ie/2f - (y')% (7) 
gAx',y') = e5iy'+ae/2)e[iw/2f-ix'n (8) 

where £^onO and are the amphtudes of the electric 
field for the TM^ q mode and the nth harmonic of the ac 
Josephson current at the surface, and the TM^ q cavity 
mode energy is degenerate for either —w/n < Xn <w/n 
or for —w/n < < w/n. The substrate factors r]j{z') 
and iriM{z') may be written as Q{z')Q{h — z') for a sam- 
ple suspended in vacuum, but for a sample on a su- 
perconducting substrate, they may be approximated as 
^gii{z')Q[K^ - (^0^][i3,[il. Using the Schelkunoff proce- 
dure of adding the contributions from these sources to the 
magnetic and electric vector potentials A(ic,t), F^{x^t) 
and/or F^ix^^t) in the radiation zone far from the mesa 
[23-25]. we have 
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where i = w^i, respectively, R = \x — x'\^ and we use the 
well-established far-field approximation 
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where r = \x\. With this approximation, the integrals in 
Eqs. ([9]) and (p!Q|) are elementary, and A{x,t) and the 
Fi{x^ t) in spherical coordinates are given in the radiation 
zone by 
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where Xn = {knW / 2) sin cos 
knW = nir/nr^ Xn a = n7r/2 - 
V = w£h, S^iO) = S^iO) = 1 
r sin 6 cos (j) -\- 6 cos 6 cos — sin ^, y = f sin 6> sin (j) + 
^ cos 6 sin (j) -\- cj) cos ^, and x^, appear in Eqs. (|5|) and 
(j6j), respectively. 

Using the Schelkunoff procedure that Ea - 

Ef = -^VxF, and = VxEp, 

Ea + Ef and H = Ha^Hf- Then, the 
differential power per unit solid angle in the radiation 
zone is given by dP/dQ = ^Re[r'^r - E x ii"*], where the 
overline denotes a time average [ii, [ill-iil. Because of 
the nearly symmetric radiation observed on both sides 
and on both ends of the mesa, with maximal radiation 
in the xz plane with 6 ^ ±30° Q , we study the incoher- 
ent mixture of the coherent radiation from both the uni- 
form ac Josephson current and the electromagnetic cavity 
sources in two models. In Model I, we preserve the von 
Neumann boundary condition that the tangential com- 
ponent of the magnetic field vanishes at the boundary, 
which becomes either Hy{x' = ±il'/2) = for standing 
waves across the width or Hx{y' = ±^/2) = for stand- 
ing waves centered symmetrically or antisymmetrically, 
depending upon £/w^ along the length. This leads to the 
respective solutions Xn = O^w/n (or yn = O^w/n) for n 
odd, and Xn = ±^^^/2n (or yn = :^w/2n) for n even, and 
we average dP/dQ over these two solutions. In Model 
II, we relax the boundary condition, and average dP/dCl 
over the allowed range of values (from —w/n to w/n) 
of either Xn or yn. This model corresponds to traveling 
waves, as the positions of the wave nodes during the time 
of the measurement can be considered to have moved an 
arbitrary amount. Then, it is straightforward to obtain 



the time-averaged dP- /dVt in the radiation zone [17|, [l 
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where j = I, II, i = wj, an{0) = \EonoS^ {0)1^(2 Zo)^ 
Zq is the vacuum impedance, kj = ccjj, and the ex- 
pressions for the C;^ ^ , I)^ ^ , and Ef^ ■ are given in the 
Appendix. For a superconducting substrate, we take 
Sl^{e)(xS^{0)occose 0,|l8|. 

The expressions listed in the Appendix exhibit the 
symmetries that the pairs D"^ ■ and and the pairs 

E^ ■ and E^ for i i' are related by interchanging Xn 
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-L cav /Ptot 

I 0.4307 0.1335 0.1235 5.431 0.491 

II 0.4241 0.2592 0.1238 3.343 0.539 

I 0.5606 0.1436 0.1250 3.320 0.328 

II 0.5567 0.2873 0.1234 1.764 0.343 



TABLE I: Least-squares fitting parameters A[ ^ and ■ for 
tiie output dPl./dQ = ^{./i(6',0) + i3(.^{.(6',0) and tiie 
standard deviations ■ obtained by assuming the waves form 
across the width and along the length in Models I and II. 

The integrated power Ptot = A{ ^fi -\- B{ ■§( •, where fi = 
1.744 for £/w = 400/77.4 and = 4.2. The portion of the 
integrated output power arising from the respective cavity 

mode IS -t cav 




FIG. 1: (color online) Polar (0 in degrees) plots of the four 
least-squares fits to the data taken in individual planes, (a) 
the xz plane (0 = 0°) normal to the mesa length, (b) The 
yz plane (0 = 90°) normal to the mesa width. Solid orange 
curves: w, I. Dashed black curves: w, II. Short-dashed olive 
curves: i, I. Dash-dotted red curves: i, II. Data are from 
Kadowaki et a/.^]. See text. 



with Yn and Xn^a with Fn,cr- The relations between the 
Model I quantities are broken into even and odd values 
of n. We note that for square samples, i/w = 1, these 
symmetries merely reflect a rotation by 90° about the z 
axis. In this case, radiation from either wave formation 
direction is equally allowed [20|. 



III. FITS TO THE AVAILABLE DATA ON A 
RECTANGULAR MESA 

We have fit the angular dependence of the radiation 
from the rectangular mesa of Kadowaki et al. [3] , using 
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FIG. 2: (color online) (a) Three-dimensional plot of the best 
fit of dPl^yj (^,0)/ dO. to the angular data of Kadowaki et al. [J] , 
obtained with the standing wave across the width, (b) Three- 
dimensional plot of the best fit of dPl^^{0, (j))/dVL to the same 
data obtained with the traveling wave across the length. 



Models I and II for the n = 1 mode, for both waves 
forming across the width and along the length, including 
the uniform ac Josephson current in all cases. We note 
that the rectangular axes were defined by Kadowaki et 
al. to be in accordance with antenna theory [23|, , with 
the X, and z axes respectively along the height, length, 
and width of mesa, so that was the angle of a vector 
from the x axis and (j) was the angle of the projection 
of the vector in the yz plane from the y axis [4]. In 
this work, we assume the standard definitions of spherical 
coordinates that the x, z axes correspond respectively 
to the directions across the width, length, and height of 
the mesa, so that is the angle a vector makes with the 
z axis, and (j) is the angle the projection of the vector 
in the xy plane makes with the x axis. There are four 
values for the fitting parameters, which we denote A\ ^ 



and ^, where i = w^£^ and j = I, 11. These A{ ^ and 
B{ ■ are fitting parameters of the theoretical functions to 
the data, 
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where the Xn{0,(j)) are given by Eq. (p!3|) and the 
Ci^.{0,(l)), D{ .{0,(I)) and E{^.{0,(I)) are given in the Ap- 
pendix. 

The n = 1 mode of a rectangular mesa with mean 
£/w = 400/77.4 ^ 5.17 is the only sample for which 
the angle dependence of the radiation in both the xz 
and yz planes is available [4]. The least-squares fits 
include the 85 data points in the xz plane (0 = 0°) 
and the 54 data points in the yz plane {(j) = 90°) with 
—90° < < 90°. The nearly symmetric distribution 
about ^ = 0° implies a rather incoherent mixture of 
the two radiation sources P, 0, Elj. In Table I, a 
list of the n = 1 fitting parameters A{ ^ and ^ for 
j = I, II and i = w,£ is given. In addition, numeri- 
cal values for the integrated components gl ^ to the out- 
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put power P{. = J^^^ sinOdO J^"" d(l)dPi .{0,(t))/dn from 
the respective cavity modes with £/w = 400/77.44 and 
rir =4.2 are given, and from the hsted Ai ^ and B{ ^ val- 
ues and the analogous integral fi = 1.744, the fraction 

Pcav/Ptot = [1 + AlJ^/{Bl^gy]-^ of the integrated 
output power arising from the respective cavity mode 
is also given. We remark that the standard deviations 
of all of the fits are very similar, so that either cavity 
wave formation direction is equally probable. However, 
the nearly identical overall best fits are for the standing 
wave across the width (Model I) and the traveling wave 
along the length (Model II), as seen in Table I. 

In Fig. 1, polar {0) plots of the best two-parameter 
fits to the experimental data of Kadowaki et at. in the 
xz and yz planes for each of the four models are shown. 
While all four models give nearly identically good fits 
to the data in the xz plane, as shown in Fig. 1(a), the 
third model with the wave across the length in Model 
I leads to a distinctly worse fit to the data in the yz 
plane, as shown by the short-dashed olive curve in Fig. 
1(b). In addition. Model II for the traveling wave across 
the width does not provide quite as good a fit in the 
yz plane as do Model I for the standing wave across the 
width and Model II for the traveling wave across the 
length. Moreover, there are some physical differences in 
the nature of the fits due to the dual-source mechanism. 
In the cases of an excited cavity mode which is either a 
standing or traveling wave across the width, nearly equal 
proportions of the radiation arise from the excited cavity 
mode and from the uniform ac Josephson current source 
[4]. On the other hand, if the excited cavity mode were 
to result in either a standing or traveling wave along the 
length of the mesa, only about one-third of the radiation 
would arise from the excited cavity mode, and about two- 
thirds of it would arise from the uniform ac Josephson 
current source, as indicated in Table 1. 

Three-dimensional plots of the predicted angular de- 
pendence of the output radiation obtained from the two 
best fits are shown in Fig. 2. As in Kadowaki et al. 
01, Fig. 2(a) is for the dual-source mechanism with the 
standing wave across the mesa width in Model I, pre- 
serving the von Neumann boundary conditions. The an- 
gular dependence of the output radiation obtained for 
the traveling wave along the mesa length in Model II is 
shown in Fig. 2(b). In this case, the wave can move 
back and forth along the length without regard to the 
von Neumann boundary conditions. As seen from Figs. 
1 and 2, it is nearly impossible to decide visually which 
fit is better, as the three-dimensional patterns formed in 
the radiation zone are nearly identical, and the standard 
deviations only differ by 7 x 10~^, which is completely 
negligible. 




FIG. 3: (color online) Spherical plots of the electromagnetic 
cavity radiation from a rectangular mesa with ijw — ^ and 
Ur — 4.2 suspended in vacuum, (a) The TM2,o raode (b) The 
TMg^g mode. 




FIG. 4: (color online) Spherical plots of the electromagnetic 
cavity radiation from a rectangular mesa with l/w — A and 
Tir = 4.2 suspended in vacuum, (a) The TMs^i mode (b) The 
TM|^2 mode. 

IV. SPECIAL CASES 

In order to guide future experiments to help determine 
which cavity mode excitation actually radiates, we now 
examine the radiation from the cavity source in some 
special situations. We are particularly interested in the 
azimuthal angles (j) = 0, 90°, for which the radiation is in 
the xz and yz planes, respectively. Let us consider the 
special case 



where q is any integer. For this special case, it is easy 
to see from Eqs. ([I7j) and ^ that for ^ = 0°, Fn = 0, 
Fn,cr = ^TT, and hence for arbitrary yn^ i ~ £ ~ ^• 
Thus, for either a pure cavity radiation source or a coher- 




FIG. 5: (color online) Spherical plots of the electromagnetic 
cavity radiation from a square mesa with i/w = 1 and rir = 
4.2 suspended in vacuum, (a) The TMf ,o mode (b) The TMf ,i 
mode. 
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ent or incoherent mixture of the two coherent radiation 
sources, for ni/w = 2q^ one would not expect to obtain 
any radiation in the xz plane from an excited cavity mode 
leading to a wave along the length of the mesa. Hence, 
for these special cases, the radiation in the xz plane is 
expected to arise from the uniform ac Josephson current 
source and/or from the excited cavity mode forming a 
wave across the width. 

There are two different examples of the special case. If 
n = 1 or some odd integer, for the special case to be re- 
alized, we must have i/w = 2qm. order that no radiation 
in the xz plane from cavity modes forming waves along 
the length could radiate. In the case n = 2 or any higher 
even value, this special case would apply to any integral 
value oi ^/w. 

On the other hand, symmetry suggests we also examine 
the yz plane, for which (j) = 90°. In this case with n£/w = 
2q and = 90°, Xn = 0. Although for n odd, Xn,a = 
n7r/2 is an odd multiple of 7r/2, for even n, Xn,a = n7r/2 
is an integral multiple of tt. As seen from Eqs. (p!5]) and 
(fTG]) . n even and ^ 90° implies M^^ = M^^ = 
for any value of x^, and the radiation from the excited 
cavity mode forming a wave across the width vanishes in 
the yz plane. Thus, examination of rectangular mesas 
with integral £/w that radiate in the TM^ q mode would 
be especially illuminating to separate out the possible 
standing wave contributions to the radiation. 

We note that when n£/w = p is din integer, /n,o = /o,p, 
so the energies of the modes with n half-integral wave- 
lengths across the width are degenerate with those of 
p half-integral wavelengths along the length. Thus, si- 
multaneous radiation from both cavity modes is possible. 
The amplitudes of the electric field for these cases may 
be written as £^onO and ^oop, which are not necessarily 
identical. In this case, the combined radiation at the nth 
mode frequency /n,o may be written as a modified version 
of Eq. (jSOl), as 
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(25) 
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given frequency: that due to the uniform ac Josephson 
current source, and that due to the excitations of the 
degenerate cavity TM^ q and TMg^ modes. Hence, a 
three parameter fit is required to determine the relative 
contribution of each source. 

We now consider the case that the rectangular cavity is 
excited into a non-degenerate TM^^ mode with m,p 7^ 
by resonance with the fundamental ac Josephson current 



angular frequency ujj 



given by Eqs. (2) 



and (3). The magnetic surface current density for the 
TM^^ mode may be written as 
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where fa{x'^y') and g^^x' ^y') are given by Eqs. ([7]) and 
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Then, the x and y components of the magnetic surface 
current density can be easily evaluated. Aside from the 
usual overall constants, in the radiation zone, the contri- 
bution to the electric vector potential from the (m,p)th 
mode is 
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xS^^{e){xM-^^ + yMy,^), (34) 
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\Eoop?/\EouO?, and the C; 
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and 



where Rn 

E^ ■ for m = n^p are given in the Appendix by setting 
kn = kno = kop. Thus, when n£/w = p, in principle 
there will be three contributions to the radiation at a 



We first preserve the von Neumann boundary condition, 
so that we average dP/dVt in the Model I procedure 



over both Xm and yp\ {dP{xm) / dQ)\ 
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[dP{{))/dQ. 



dP{w/m)/dQ] for m odd and {dP{xm) / dQ)i = 
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FIG. 6: (color online) Spherical plots of the electromagnetic 
cavity radiation from a square mesa with l/w = 1 and rir = 
4.2 suspended in vacuum, (a) The TM|,o mode (b) The TM2,2 
mode. 



\[dP{-w/2m)/dQ. + dP{w/2m)/dQ\ for m even. Sim- 
ilarly, {dP{yp)/dn)i = ^[dP{0)/dQ + dP{i/p)/dQ] 
for p odd and {dP{yp)/dn)i = \[dP{-i/2p) / dVt + 
dP{£/2p)/drt] for p even. In Model II, we average over 
all allowed values of —w/m < Xm < w/m and —i/p< 
Vp ^ ^/p^ and the results are given in the Appendix. 

Three-dimensional plots of the radiation from the 
TM2 mode and the degenerate TMg g mode for i/w = 4 
and nr = 4.2 for a sample suspended in vacuum as calcu- 
lated in Model I are shown in Fig. 3. Spherical plots of 
the TM| 1 and TM| 2 modes with £/w = 4 and = 4.2 
are shown in Fig. 4. We note that the radiation from the 
TM2 mode has a maximum in the xz plane, whereas 
the radiation from the degenerate TMg g mode vanishes 
in the xz plane and has a maximum at a finite value 
in the yz plane. For the TM^ i and TM^ 2 modes with 
slightly higher energies, the radiation patterns shown in 
Fig. 4 indicate that there is also no radiation in either 
the xz and yz planes. 

In addition, in anticipation that angular measurements 
of the radiation from square mesas might soon be made 
[20] , we plotted the anticipated radiation from the cavity 
excitation of a square mesa with i/w = 1 and = 4.2 
using Model I for a sample suspended in vacuum. In 
Fig. 5(a), we plotted the radiation for the fundamental 
TMf mode. This mode has maximal radiation in the 
xz plane. The degenerate TMg ^ mode exhibits radiation 
that is precisely the same as in Fig. 5(a) but rotated by 
90° about the z axis. In Fig. 5(b), the radiation from 
the cavity TMf ^ mode is shown. The radiation from this 
mode is invariant under rotations of 90° about the z axis, 
as expected from symmetry arguments, is small in both 
the xz and yz planes, and vanishes at ^ = 0°. In Fig. 6, 
the radiation of the TM^ and TM^ 2 modes of a square 
mesa with rir = 4.2 are shown. As seen from Fig. 6(a), 
the radiation from the cavity TM^ q mode vanishes in the 
yz plane, and symmetry thus requires that the radiation 
from the cavity TMg 2 mode vanishes in the xz plane. 
The cavity radiation from the TM^ 2 mode vanishes in 
both the xz and yz planes, as shown in Fig. 6(b). 

To illustrate the effect of a combination of the radia- 
tion from the degenerate cavity modes, in Fig. 7, three- 
dimensional plots of the combined cavity radiation from 




-JUT 




FIG. 7: (color online) Spherical plots of the combined cavity 
radiation from the degenerate TM2,o and TMq^s cavity modes 
in rectangular mesas with integral l/w — A and Ur — 4.2. 
(a) 1^020 1/1^008 1 = 1 for a sample suspended in vacuum, (b) 
1^020 1/1^008 1=0.1 for a sample on a superconducting sub- 
strate. 




FIG. 8: (color online) Spherical plot of the combined cavity 
radiation from the degenerate TM^ q ^"^^ TMo,n cavity modes 
in square mesas with l/w = 1 and rir = 4.2 suspended in 
vacuum, (a) The combined TMf q and TMg 1 modes with 
l^oiol/l^ooil = 1. (b) The combined TM^^q and TMg,2 modes 

with 1^020 1/1^002 I = 1. 



the TM2 and TMg g modes of a sample with l/w = A 

and Ur = 4.2 are shown. In Fig. 7(a), |£^oo8|/|^02o| = 1 
for a sample suspended in vacuum. For an equal weight 
of the two contributions, the largest radiation occurs 
for a narrow range of values in the yz plane. In 
Fig. 7(b), the analogous combined cavity radiation with 
1^008 1 / 1^^020 1 = 0.1 for a sample on a superconduct- 




FIG. 9: (color online) Spherical plots of the uniform ac 
Josephson current radiation from rectangular mesas with in- 
tegral i/w and Ur = 4.2 suspended in vacuum, (a) i/w = 4 
with the frequency locked onto that of the TM2,o cavity mode, 
(b) i/w = 1 with the frequency locked onto that of the TMf 
cavity mode. 
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ing substrate is pictured. In Fig. 8, three-dimensional 
plots of the equally combined radiation from the TM^ q 
and TMq ^ modes with = 4.2 for a square sample 
suspended in vacuum are pictured. In Fig. 8(a), the 
combined TMf q and TMg ^ cavity mode radiation with 

1^010 1 / 1^001 1 = 1 is shown. In this case, there is no 
azimuthal anisotropy, and the axial anisotropy is rather 
weak. The radiation is a maximum at ^ = 0°, as for cylin- 
drical cavities in the cylindrical TMf ^ mode lH, 20 1. In 
Fig. 8(b), the combined radiation from the cavity TM^ q 

and TMq 2 modes with |^oo2|/|^02o| = 1 for the square 
sample suspended in vacuum is also shown. For these 
modes, the radiation exhibits C^y point group symmetry, 
but is highly anisotropic, not only in the axial direction, 
but also azimuthally. 

However, as it has been established that the radiation 
from existing mesas of Bi-2212 arises from a dual-source 
mechanism, with the primary radiation source being the 
uniform ac Josephson current [4, 20], and further sup- 
ported by the four fits to the existing data shown in 
Figs. 1 and 2, it is important to remember that the 
predicted electromagnetic cavity radiation displayed in 
Figs. 3 - 6 for rectangular mesas with £/w = 4^1 is only 
part of the predicted radiation. In addition, there will 
be the radiation from the uniform ac Josephson current 
source with the same frequency. In Fig. 9, we displayed 
the radiation patterns corresponding to those patterns in 
Figs. 3-6. In Fig. 9(a), the radiation from the uniform 
ac Josephson current source in a rectangular mesa with 
£/w = 4 and rir = 4.2 suspended in vacuum with the 
frequency locked onto that of the cavity TM^ q mode is 
shown. This pattern is exactly the same as that for the 
frequency locked onto the TMg g mode, as those frequen- 
cies are identical. Moreover, the patterns for the slightly 
different frequencies locked onto the TM^ i and TM^ 2 
modes are almost indistinguishable from that shown in 
Fig. 9(a), and hence are not shown. Thus, although the 
radiation patterns from the excitation of TM^^ electro- 
magnetic cavity modes with p = 0, 1, 2 are extremely dif- 
ferent, the patterns predicted to arise from the primary 
radiation source, the uniform ac Josephson current, are 
nearly identical. We note that for both sources, the ra- 
diation at ^ = 0° vanishes, unlike the case of the electro- 
magnetic cavity radiation for the rectangular mesa with 
mean i/w = 400/77.4, the fits to which are displayed in 
Figs. 1 and 2. In addition, we note that if the sample 
were placed upon a superconducting substrate, the ra- 
diation patterns would all vanish at 6> = 90°, as in Fig. 



on patte 
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7(b) 

In Fig. 9(b), we show the predicted radiation from the 
uniform ac Josephson current source from a square mesa 
{i/w = 1) with Ur = 4.2 suspended in vacuum, with the 
frequency locked onto that of the TMf q and the degen- 
erate TMq I mode. Unlike the highly asymmetric corre- 
sponding electromagnetic cavity radiation shown in Fig. 



5(a) which has a maximum at ^ = 0° and is strongly ^- 
dependent, this radiation is independent of the azimuthal 
angle (/), and vanishes at ^ = 0°, just as for the ac Joseph- 
son current source radiation, the frequency of which is 
locked onto that of the cylindrical cavity TMf ^ mode 
(2q| . The lack of azimuthal anisotropy was also present 
for the combined radiation of the cavity TMf q and TMg ^ 
modes with |^ooi |/|^oio| = 1 pictured in Fig. 7(a), but 
in that case, the radiation was a maximum at ^ = 0°, 
and in this case, it vanishes at 6> = 0°. Both cases thus 
are essentially the same as for the corresponding cases 
for the uniform ac Josephson radiation locked onto the 
cylindrical cavity TMf ^ mode [20] . For the frequency 
locked onto that of the rectangular TMf ^ mode, the ra- 
diation emitted from the uniform ac Josephson current 
source is nearly identical to that pictured in Fig. 9(b), 
so it is not shown for brevity. In addition, the radia- 
tion from the ac Josephson current source the frequency 
of which is locked onto either that of the cavity TM^ q 
mode or the cavity TM^ 2 mode is also independent of 
azimuthal angle ^, so those figures are also not shown 
for brevity. At those frequencies, both radiation sources 
lead to a vanishing output at ^ = 0°, but each cavity 
source radiation is highly ^-dependent, unlike the uni- 
form ac Josephson current source radiation. As for the 
other rectangular mesas studied, if the square mesa were 
placed atop a superconducting substrate, both radiation 
sources would lead to a vanishing output at ^ = 90°. 



V. DISCUSSION AND SUMMARY 

Since the publication of the Wang et al. paper [2l|, 
there has been an ongoing debate as to whether the ra- 
diation emitted from rectangular mesas originates partly 
from the standing wave across the width of the mesa or 
partly from excited waves of the same wavelength \ = 2w 
that appear to form along the length of the rectangular 
mesa. Since experiments on a large number of rectan- 
gular mesas have determined that the wavelength of the 
radiation \ = 2w [l-ll], it is natural to suppose that the 
cavity mode to which the ac Josephson current couples is 
the TMf Q mode. However, laser measurements have sug- 
gested that the standing waves forming along the length 
of the mesa appear to have the same wavelength A ^ 2k;. 
Since the mesas in which these standing waves have been 
detected have since been shown to emit radiation, it is 
important for our understanding of the radiation mecha- 
nism to determine which, if either, of these possibilities is 
correct. We therefore re-examined the published angular 
dependence of the radiation emitted from a rectangu- 
lar mesa [1]. Although one might suppose that it would 
be easy to distinguish these models from the calculated 
output radiation, for the dual-source mechanism with an 
arbitrary mix of radiation from the uniform ac Josephson 
current source and the radiation from the excited electro- 
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magnetic cavity source, we found that both scenarios are 
equahy probable: it is impossible from the existing data 
to distinguish whether the cavity source forms with the 
wave across the width or along the length of the rectan- 
gular mesa. This conclusion is based upon two models of 
the coherent cavity and ac Josephson radiation sources 
contributing incoherently to the output, as required for 
any reasonable fit to experiment, since the data for pos- 
itive and negative values were nearly symmetric both 
in the xz and yz planes [J]. The equally probable over- 
all best fits to the data of Kadowaki et al. from a mesa 
with mean l/w = 400/77.4 are for the standing wave 
forming across the mesa width in Model I and for the 
traveling wave forming along the sample length in Model 
11. The first case is consistent with the analysis of Kad- 
owaki et a/.Q, preserving the von Neumann boundary 
conditions. In the second case, the traveling wave is free 
to move along the length of the sample with no preferred 
nodal position. 

In addition, there is a physical difference in the relative 
importance of the cavity portion of the overall output 
radiation of the two radiation sources. If the excited 
cavity waves were to form across the sample width, about 
half of the total radiation would arise from each source. 
But if the excited cavity waves were to form along the 
sample length, about two-thirds of the radiation would 
arise from the ac Josephson current source. 

We remark that there are many ways to do the sym- 
metric averaging of the cavity mode radiation in the 
presence of the radiation from the uniform ac Joseph- 
son current source. In this work, we only used Models 
I and II. For Model I with the wave along the length, 
we assumed the wave was centered either symmetri- 
cally, for l/w slightly larger than an even integer, or 
asymmetrically, for l/w slightly larger than an odd in- 
teger, as for the sample under study. Then, we assumed 
{dPi^e{yn)/dn)i = l[dP{0)/dQ^dP{w/n)/dQ] for n odd 
and ^[dPi/{—w/2n)/dft^dPi^£{w/2n)/dQ] for n even. 
This results in a contribution to the output from the 
two-source mechanism that is symmetric about ^ = 0°, 
as in experiment [3]. It also preserves the von Neumann 
boundary condition. We note, however, that this Model 

I fit for a standing wave along the length led to the worst 
fit of the four models, but the Model II fit for a traveling 
wave along the length led to a tie for the best overall 
fit to the data. Hence, some other averaging of yi that 
preserves the even symmetry of the combined output 
of the two-source mechanism could lead to an improved 
overall best fit. For example, one could in principle treat 
the combined equal outputs from yi and yi -\- as in 
the experiment, treating yi as a fitting parameter in ad- 
dition to Al/ and Si/, and perform a three-parameter fit 
to the data. Such a three-parameter fit would likely lead 
to a lower standard deviation than that for the Model 

II fit, in which yi was averaged over all allowed values. 
This optimum fixed value of yi would lead to a standing 



wave along the mesa length. However, only about 1/3 
of the combined output arose from the excited cavity 
mode in the Models I and II fits. Hence, modifying the 
fit with a third adjustable parameter that improves the 
minor output fraction without making some comparable 
adjustment to the major (uniform ac Josephson current 
source) output fraction could be unreliable. A more re- 
liable three-parameter fit might therefore arise from the 
inclusion of a non-uniform correction to the ac Josephson 
current source jl8|] . Hence the question as to whether the 
optimum fit for a wave along the rectangle length would 
be standing or traveling is unresolved. 

We then considered the possibility of investigating 
mesas with different £/w ratios. We found that for 
ni/w = 2g, where q is any integer and n is the index 
of the TM^ Q rectangular cavity mode, the electromag- 
netic cavity radiation from a standing wave along the 
length of the mesa (in the degenerate TMg 2^ mode) van- 
ishes in the xz plane. In those cases, any observed ra- 
diation in the xz plane would have to arise from either 
the TM^ Q electromagnetic cavity mode or the uniform 
ac Josephson radiation. In addition, for n even, the ra- 
diation from the cavity TM^ q mode with a wave across 
the width vanishes in the yz plane. We presented figures 
of the predicted radiation for a rectangular sample with 
i/w = 4 and for a square sample, both suspended in vac- 
uum, to illustrate these points. Since for integral i/w, 
the radiation can occur from either or both of the two 
degenerate cavity modes, we pictured examples of the 
combined cavity radiation from both modes, which has 
reduced azimuthal anisotropy, especially when i/w = 1. 
We strongly encourage further experiments on samples 
with even-integral i/w ratios to confirm these predic- 
tions. It would be very interesting to compare samples 
with even-integral £/w values with those having slightly 
different £/w values, to see if the radiation pattern is 
substantially different. 
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rials Nanoarchitectonics) project (NIMS) and Strategic 
Initiative category (A) at the University of Tsukuba. 



APPENDIX 

First, we list the results of our calculations for the ra- 
diation from a general cavity mode for an arbitrary i/w 
ratio. For Model I, we have 

Dl. = Bl„ (38) 

Eli = 2An,iBn,i (39) 
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for i = w^i, where 
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Then, we hst the results for the radiation from the 
TM^^ mode for a rectangle of length ^ and width w. 
Then, in Model I, we have 



2 ' -mp ^ 
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In Model II, we find 

clip = 2sin2(f 



(53) 



(54) 
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■)• (55) 



Then, the differential power per unit solid angle emitted 
from the TM^^ cavity mode is proportional to 



^ |5^,Wp(c^, + I^^^,-sin^^[C^^cos^ 



sin (j) - El^^ sin (j) cos ( 



(56) 



for j = I, II. 
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